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Spontaneously em erg ent chirality is an issue of fundamental importance across 
the natural sciences (l(|. It has been argued that a unidirectional (chiral) rota¬ 
tion of a mechanical ratchet is forbidden in thermal equilibrium, but becomes 

.. . n 

possible in systems out of equilibrium p||. Here we report our finding that a 
topologically nontrivial spin texture known as a skyrmion - a particle-like ob¬ 
ject in which spins point in all directions to wrap a sphere (3[| - constitutes 
such a ratchet. By means of Lorentz transmission electron microscopy (Lorentz 
TEM) we show that micron-sized crystals of skyrmions in thin films of Cu 20 Se 03 
and MnSi display a unidirectional rotation motion. Our numerical simulations 
based on a stochastic Landau-Lifshitz-Gilbert (LLG) equation suggest that this 
rotation is driven solely by thermal fiuctuations in the presence of a tempera¬ 
ture gradient, whereas in thermal equilibrium it is forbidden by the Bohr-van 
Leeuwen theorem jj, 1^. We show that the rotational fiow of magnons driven 
by the effective magnetic field of skyrmions gives rise to the skyrmion rotation, 
therefore suggesting that magnons can be used to control the motion of these 
spin textures. 

The formation of triangular arrays of skyrmions in chiral-lattice magnets under an allied 


magnetic fie 

Fei_3;Co3;Si 


d was theoretically predicted (6|-l8| a^ 
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13|, FeGe [l^, and Cu20Se03 


15 


experimentally observed in MnSi [9, 10|, 


17|. The magnetization in such crystals 


is antiparallel to the magnetic field B at the center of each skyrmion and is parallel to B 
at its periphery (see Figs.la-c). Recently, real-space images of nanosized skyrmions have 
been successfully obtained by the Lorentz TEM. While scanning temperatures and magnetic 
fields in the experiments, we have encountered a peculiar dynamical phenomenon, that is, 
in a wide temperature interval excluding only the lowest temperatures, micro-scale regions 
of the skyrmion crystal (SkX) show, in addition to Brownian motion, a clearly discernible 
unidirectional rotation. 

In Supplementary Movies 1 and 2, we show examples of this phenomenon observed in thin- 
plate (~50nm-thick) specimens of MnSi and Cu20Se03, respectively. These compounds are 
chiral-lattice magnets with a common space group P2i3. In zero field they undergo a phase 
transition from paramagnetic to helimagnetic phase at 29.5 K and 60 K, respectively. The 
helix period \m is, respectively, 18 nm and 50 nm. Thin-plate (<100 nm-thick) specimens 
of these compounds host stable SkX phases with the skyrmion lattice spacing ~ ;^Am- The 
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field strength for the Movies 1 and 2 is 175 niT and 65 niT, respectively, and B is applied 
in the negative 2 ;-direction perpendicular to the plate plane: B || —z. The observed (static) 
Lorentz-TEM image of SkX in MnSi specimen is displayed in Figs, dt and d. In Figs. [I^-l, we 
exemplify time evolution (snapshots) of Fourier components of the magnetic configurations 
where a hexagon composed of six Bragg peaks rotates clockwise, indicating the clockwise 
rotation of the SkX domains. 

The rotation rate depends on the irradiation density of the electron beam of the micro¬ 
scope (Supplementary Information). The skyrmion rotation is observed only above a critical 
irradiation density and the rotation rate grows as the density increases. This indicates that 
this rotation is a non-equilibirum phenomenon induced by the electron beam. The similarity 
between the chiral rotations of SkX micro-domains and their hexagon-shaped Fourier com¬ 
ponents observed in MnSi (Supplementary Movie 1) and Cu 20 Se 03 (Supplementary Movie 
2) is remarkable in view of different origins of magnetism in these two compounds (MnSi is a 
metal, while Cu 20 Se 03 is an insulator) and the difference in skyrmion parameters, e.g. the 
transition temperature and the SkX lattice spacing. Hence, such a ratchet motion should 
be viewed as a generic feature of skyrmion systems. Then the question arises why and how 
it occurs. 

A circular magnetic held induce by the electron beam of the Lorentz TEM is estimated 
to be hve orders of magnitude smaller than the geomagnetic held, and thus cannot cause 
the rotation. Recently it was demonstrated that a spin-polarized electric current parallel 
to the sample plane can drive similar rotations of SkX domains in the presence of thermal 
gradient 20, l2]|. In our case, however, a possible ehect of electric currents can be excluded 
since the electron beam of Lorentz TEM is three orders of magnitude smaller than the 


threshold current of 10^-10® A/m^ for the current-driven skyrmion motions 


3 - 23 | 


and, in 


addition, its direction is perpendicular rather than parallel to the him, which further reduces 
spin-torque ehects. Furthermore, the energy of the electrons is very high and the interaction 
with the spins in the specimen is very weak except through the magnetic held induced by 
the magnetization. This conclusion is corroborated by the fact that the skyrmion rotation 
is also observed in the insulating Cu 20 Se 03 . Thus, the unidirectional rotation is apparently 
induced by thermal ehects. 

In order to clarify the nature of this phenomenon, we performed numerical simula¬ 
tions. The SkX phase in chiral-lattice magnets is described by a classical Heisenberg 
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model on the two-dimensional square lattice 18|, which contains ferromagnetic-exchange 


and Dzyaloshinskii-Moriya 


DM) interactions as well as the Zeeman coupling to B={0, 0, 


Bz) normal to the plane 19|]. The Hamiltonian is given by, 


T-i = —J rrii ■ rrij — D x rrii^s: • x -|- x rrii^^ ■ y) 

<i,j> i 

-BE 


m. 


( 1 ) 


where the magnetization vector is dehned as mi = —Si/h with the norm m=|mj|. The 
spin turn angle 6 in the helical phase is determined by the ratio D/J and for D/J=0.27 
used here, 9=11° corresponding to the period of ~33 sites. With a typical lattice constant 
of 5 A, this gives Am ~17 nm, which is comparable to the helical period ~18 nm in MnSi 
(our conclusions, however, are not affected by the choice of the value of D/J). Our Monte- 
Carlo analysis shows that the SkX phase emerges in the range of 1.688 x 10“^ < Jm < 
5.67 X 10“^ between the helical and ferromagnetic phases. 

In this study, we treat SkX confined in a micro-scale circular disk with the diameter 
2R = 137 sites, as shown in Fig. [2^. We impose the open boundary condition and sim¬ 


ulate thermally-induced c 
stochastic LLG equation 

dmi 1 


yiiamics of this skyrmion microcrystal by numerically solving a 


24 l25| using the Heun scheme 


26l |. The equation is given by 


dt 1 -|- 


G 


mi X {Bf + + —mi X [mi x {Bf + ^f{t))] 

m 


( 2 ) 


where gg is the Gilbert-damping coefficient and Bf^ = is the deterministic held. 

The Gaussian stochastic held $f{t) describes ehects of thermally huctuating environment 
interacting with mi, which satishes = 0 and (C\/ 3 (^)'Cj,A(^)) = — s), 

where f3 and A are Cartesian indices. The huctuation-dissipation theorem gives a relation 
between k, and temperature T\ k = aQk-BT/m 26|]. The initial spin conhguration (Fig. |2^) 


is prepared by the Monte-Carlo thermalization at low temperature and by further relaxing 
it in the LLG simulation at T=0. Starting from this initial conhguration, we generate the 
random numbers corresponding to the stochastic force ^f{t) and solve equation (E]). In what 
follows we use units in which the lattice constant a = 1, the exchange energy J = 1, the 
Boltzmann constant = 1 and h = 1. 

In thermal equilibrium we only hnd Brownian motion of skyrmions and no unidirectional 
rotation. We then include a radial temperature gradient to examine whether it can give rise 
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to chiral rotation. In the Lorentz-TEM experiment, the electron beam is irradiated onto a 
thin-plate specimen, which inevitably raises temperature of the beam spot with respect to 
the outer region, resulting in the temperature gradient as shown in Fig. [2]b. We consider a 
constant temperature gradient with T =Tq at the edge to T=To + AT at the center in the 
circular-disk system, i.e., —dT/dr=AT/R. 

We display snapshots of the calculated real-space magnetization dynamics at selected 
times in Figs. |3^-d, and the trajectory of a selected skyrmion indicated by solid rectangles 
in Fig. [3li. We hnd persistent rotation of SkX (Supplementary Movies 3 and 4). Figures |3^-1 
show the time evolution of Fourier transform of the spin structure - the rotating hexagon 
composed of six Bragg peaks (Supplementary Movie 5). In the simulation, we apply B || —2 
in accord with a setup of the Lorentz-TEM experiment, and hnd clockwise rotations in 
agreement with the experimental observations. Remarkably, this unidirectional rotation is 
driven purely by the thermal gradient because no other motive forces are considered in our 
simulation. 

This nonreciprocal dynamics of SkX can be traced back to the algebra of spin operators, 
which determines the direction of spin precession in an applied B. Equations of motion for 
the center-of-mass coordinates of a skyrmion have the form (Supplementary Information) 


MY + acTY + AttQSX = 


MX + acTX - MQSY = 


dU 

dY 

au 

ax 


+ 47rQ 


( 3 ) 


where M is the skyrmion mass, Q is the topological charge {Q = —1 for < 0), 
aQ is the Gilbert-damping constant, T S.STTttS', U is the external potential, and 
jmagnon _ jmagnon^ jmagnon^ jg maguou Current density dehned in Supplementary In¬ 
formation, which drives skyrmion motion through the spin transfer torque. The extension 
to a many skyrmion problem is straightforward. 

First we note that, if one replaces the magnon current by the stochastic Langevin force. 


es in an external 
which forbids 


the dynamics of skyrmions becomes equivalent to that of classical partic 
magnetic held. In this case, the famous Bohr-van Leeuwen theorem J, 
orbital magnetism of classical particles in thermal equilibrium, precludes the spontaneous 
rotation of skyrmions, in agreement with our numerical results. 

Next we discuss how a nonzero temperature gradient can induce a persistent rotation. 
One possible scenario is that the gradient dT/dr modihes the radial distribution of skyrmions 
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resulting in a non-vanishing radial force — {dU/dr), which according to equation ([3]) gives 
rise to a nonzero angular velocity. This, however, is forbidden, as the Bohr-van Leeuwen 
theorem can be generalized to the case of local equilibrium with a spatially inhomogeneous 
temperature T{r). What is required, is a nonequilibrium state with a heat flow forming 
a heat engine, in which an amount of heat Qi transferred from the high-temperature side 
(T = Ti) is partly transformed into work W, while the remaining heat Q 2 = Qi — W is 
absorbed on the low-temperature side (T = T 2 < Ti). The ratchet rotation requires W > 0 
and the engine efficiency rj = W/Qi is less than 1 — T 2 /T 1 . 

The underlying microscopic mechanism involves the flow of the spin energy. There are 
two possible heat carriers: skyrmions and magnons. The dimensional analysis of the ratchet 
rotation frequency u due to the thermal motion of skyrmions gives u ~ i (-f) (ta ‘h* 
J = h = kB = a = i- units), which is too small to explain the results of numerical sim¬ 
ulations. Magnons, on the other hand, are much more efficient in driving the rotation of 
skyrmions through the spin transfer torque 27|]. Importantly, a skyrmion induces an effec¬ 
tive magnetic field = —m ■ (cTm x dyin) with the total flux Jdrhz = 4:nQ, which exerts 
the Lorentz force on magnons 28(] (Supplementary Information). The skew scattering of 
magnons off skyrmions gives rise to the Topological Magnon Hall Effect: in addition to the 
thermally-driven magnon current in the radial direction, (~w)’ ^^ere is 

a current direction transverse to the temperature gradient 

corresponding to the counterclockwise rotation of the magnon gas. 

Figure 0] displays the result of simulations for the magnon current density. We show a 
real-space map of the magnon current density at a selected time in Fig. 0^ where the arrows 
point in the current directions while their lengths represent the current amplitudes. In this 
map, the currents may seem to flow in random directions, but we find a net positive value 
of the time-averaged quantity x shown in Fig. 0 )d. Here 

the vector connects the center of the disk and the ith. site, and the summation goes 

over the sites with |rj| > 20. The positive sign indicates that the magnon current flows 
in the counterclockwise direction. The magnitude of the transverse magnon current in our 
numerical simulations, j™^snon ^ 10 “^^ which 


large magnon Hall conductivity, ~ 10 


or —dT /dr ~ 10 ^ corresponds to a rather 


29|. 


The skew magnon scattering off skyrmions exerts a reaction force on skyrmions in the 
negative d-direction. The reaction force appears in the right-hand side of Eqs.(0]), e.g., the 
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force Fy equals a product of the flux of the effective magnetic field AttQ and the magnon 
current From equations of motion ([3]), we obtain the estimate for the rotation rate 

of SkX (Supplementary Information) 


jy 


^magnon 

ttR 


-5 X 10 


-5 


(4) 


for 10“^. The minus sign corresponds to clockwise rotation of skyrmions. Equation 

(14|) shows that the experimentally observed clockwise rotation of skyrmions is a consequence 
of the anticlockwise rotation of magnons. The estimate (jl]) is in good quantitative agreement 
with the numerical result for the SkX rotation rate, —3 x 10“^. 

We thus showed that the magnon current induced by the temperature gradient is deflected 
by the emergent magnetic field of skyrmions, which in turn gives rise to the rotation of SkX 
through the spin-transfer torque. Since the sign of J^agnon jg governed by the sign of dT/dr, 
the rotation of SkX should be reversed upon the sign reversal of the temperature gradient, 
which is indeed what we find in our simulation (see Fig. H}:: and Supplementary Movie 6). 
Also the rotation direction becomes reversed upon the sign reversal of magnetic field Bz but 
not upon the sign reversal of the DM parameter D as seen in Fig. [Hi because the former 
changes the sign of but the latter does not. This shows that skyrmion-magnon 

interactions and thermal spin fluctuations provide a key to understanding of the observed 
chiral rotation of skyrmions. 

The proposed physics behind the observed rotation is distinct from the Skyrmion Hall 
effect discussed in a recent paper by Kong and Zang [27|. They theoretically proposed that 
a longitudinal skyrmion motion due to the magnon current along the thermal gradient is 
accompanied by a small transverse motion due to the Gilbert damping. This Skyrmion 
Hall effect necessarily requires the longitudinal motion. In our case, however, the skyrmion 
motion in the radial direction (parallel to the temperature gradient) is forbidden due to 
the geometrical conhnement. Hence, no Skyrmion Hall effect is possible and the topological 
magnon Hall effect is the only source of the observed skyrmion rotation. The reaction force 
from the magnon current deflected by the effective magnetic held of the topological skyrmion 
texture drives the peculiar chiral motion. 

Note that the time scale of the rotation is microseconds in the simulation for J~1 meV, 
while it is a few seconds in the experiment. This discrepancy can be related to the strong 
sensitivity of the rotation rate to the shape of the boundary of the SkX, the magnitude of the 
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temperature gradient and sample inhomogeneities, such as impurities and defects. Thus we 
observe no rotation in a rectangle-shaped system within a realistic simulation time (limited 
by a few milliseconds), apparently because the large friction between the SkX and the system 
edges makes the rotation rate extremely slow. In the circular-disk system, the rotation rate 
decreases as AT decreases, and vanishes when AT=0. In addition, the rotation is less 
pronounced and its rate is lower in a larger-sized disk, indicating the absence of rotation in 
the thermodynamic limit. 

To summarize, we have found experimentally and explained theoretically that micron¬ 
sized skyrmion crystals behave as the Feynman’s ratchet. In the presence of the radial 
magnon flow, skyrmions exhibit persistent rotation in the direction determined by the sign 
of the topological charge of the skyrmions. The physical origin of this unusual phenomenon 
can be traced back to the chiral nature of spin dynamics. Our hnding shows how thermal 
spin fluctuations can be harnessed to control topological spin textures by irradiating them 
with light and electrons. The manipulation of skyrmions with magnons can be used to build 
all-spin memory and logic devices with low dissipation losses by replacing the charge current 
by magnon current especially in the insulating magnets j^|. 

Methods Summary 

The single crystal samples of MnSi were grown by the floating zone technique, while 
those of Cu 20 Se 03 were grown by the chemical vapor transport method. For the real-space 
imaging of spin textures, their thin specimens with thickness of ~50 nm were prepared by 
mechanical polishing and subsequent argon-ion thinning with an acceleration voltage of 4 kV. 
All experiments were performed with a transmission electron microscope (JEM2100F, JEOL) 
at an acceleration voltage of 200 kV. Images of the SkX were obtained at the over-focused 
Lorentz-TEM mode. The movies of Lorentz-TEM image were taken with the exposure time 
of 50 ms and the frame rate of 18 fps (frame per second). A liquid Helium cooling holder 
was utilized to investigate the T dependence, by which T at the specimen can be controlled 
from 6 K to 300 K. The electron beam strength is 2.7 x 10^ Am“^ (3.6 x 10^ Am“^) for 
the MnSi (Cu 20 Se 03 ) specimen. The magnetic held of 5=175 niT (65 niT) parallel to the 
electron beam was applied, and the measurement was performed at 8 K (40 K) for MnSi 



(CuaOSeOa). 
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Figure Legends 


Fig.l I Experimentally observed Skyrmions and unidirectional rotation of 
micro-scale skyrmion-crystal domains in MnSi. a, Schematic figure of skyrmion in 
chiral-lattice magnets, b, Real-space Lorentz-TEM image of single skyrmion in MnSi. The 
colour map and arrows represent in-plane components of magnetizations, c, Real-space 
Lorentz-TEM image of the static skyrmion crystal in MnSi in which skyrmions are packed 
to form a triangular lattice, d, Over-focused Lorentz-TEM image of the skyrmion crystal in 
MnSi. e-1, Time profiles of Fourier transforms of temporally changing magnetic structures 
in the skyrmion crystal observed in MnSi by the Lorentz TEM, which show a clockwise 
rotation. Open circles and solid lines are guides for the eyes. 

Fig. 2 I Setup of the numerical simulation, a, Magnetic configuration of the skyrmion 
microcrystal confined in a circular-shaped disk at T=0 where the in-plane magnetization 
components at sites iy) are indicated by arrows when mod(i 3 ;, 2)=mod(iy, 2)=0, while 
distribution of the magnetization 2 ;-axis components, rrizi, is shown by a color map. b, 
Schematics of induced thermal gradient by the electron-beam irradiation in the Lorentz 
TEM experiment. 

Fig.3 I Simulated thermally driven rotation of the skyrmion microcrystal. 
a-d, Snapshots of simulated temporally changing distribution of the magnetization z-axis 
components rrizi at t=1.9xl0^ (a), t=2.6xl0^ (b), ^=3.3x10“^ (c), and t=4.0xl0^ (d), 
which show a clockwise rotation. Here the time unit is h/J. Also shown in d is the 
trajectory of a selected skyrmion indicated by rectangles, e-1, Time evolution of Fourier 
transforms of the magnetic structure in the reciprocal space —0.27r < k^ < 0.27r and 
—0.27r < ky < 0.27r for m=l, B^/J = —0.03, aG=0.01, /cbTo/J= 0.1, and /cbAT/J= 0.006 
at selected times between t=1.90xl0^ and t=4.35xl0^ with constant time intervals of 
At=0.35xl0^, which also show clockwise rotation. 

Fig. 4 I Simulation of the magnon current density, a, Spatial distribution of 
the magnon current density ^ selected time where lengths of the arrows are 

proportional to the local amplitudes of j™a,gnon Time evolution of the quantity 
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^ which shows a finite and positive component of in the 

circnmferential direction. Here is the vector connecting the disk center and the ith site, 
and the snmmation runs over the sites with |rj| > 20. The sum of the absolute values, 
that is, Yl'i ^ is plotted to show that the deviation of the above quantity 

from zero is meaningful. The calculation is performed for m=l, J = —0.03, aG=0.01, 
A:bTo/J=0.1, and /cbAT/ J=0.006. c, Simulated number of rotations 6/2ti as functions of 
time in the cases of thermal equilibrium (AT = 0), positive T gradient (AT > 0), and neg¬ 
ative T gradient (AT < 0). Here the negative (positive) slope of the plot means a clockwise 
(counterclockwise) rotation of the SkX. The result shows that the rotation sense is reversed 
upon the sign reversal of AT. d, Simulated 9/2t[ as functions of time for four kinds of 
combinations of signs of and the DM parameter D for /cbAT/J=-|- 0.006. The rotation 
sense changes depending on the sign of Bz but not on the sign of D. 
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